
  x = [0  .906281812634947  -.640769438625687  .640806205898312  .000105490680686318  -.906268698302751  .640846171806826  -2.4376976446418e-005  -.640913478668817  -6.34307633893959e-005  ];

  y = [0  0.  .640881974571767  -.640875894070397  .906286487817418  .00013131435482042  .640758959816975  -.906264855773148  -.640815543815499  -.000108398405231571  ];


  z = [1  .422673959556436  -.422711274030294  -.422664754719007  .422663921894554  .422702057285244  -.422781432296648  .422710315223367  -.422593601069731  -.999999992113162  ];

  plot3(x,y,z,'rs','Markersize',15);
  sphere;
  axis equal;
  hold on;
  plot3(x,y,z,'rs','Markersize',15,'Markerfacecolor','black');
  hold off;
  
  s = 0;
  
  for i = 1:length(x)
      for j = i+1:length(x)
          s = s + (x(i)-x(j))^2 + (y(i)-y(j))^2 + (z(i)-z(j))^2;
      end
  end
  s